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Abstract:  In  this  paper,  a  position  regulation  control  strat¬ 
egy  is  developed  for  a  magnetic  levitation  system  operating  in 
the  presence  of  a  bounded,  nonlinear ,  periodic  disturbance.  The 
proposed  controller  utilizes  a  saturated  control  force  input  in  con¬ 
junction  with  a  learning  based  disturbance  estimator  to  asymp¬ 
totically  regulate  the  target  mass  to  a  desired  set  point  position 
despite  the  actuator’s  unidirectional  limitation  of  exerting  only 
an  attractive  force  on  the  target  mass  (i.e.,  the  control  input  can 
only  generate  an  attractive  force  while  the  earth ’s  gravitational 
field  is  utilized  to  produce  the  repulsive  action).  In  addition, 
the  control  development  only  requires  the  disturbance  input  to  be 
bounded  and  periodic  in  nature.  Simulation  results  are  included 
to  illustrate  the  performance  of  the  control  strategy. 

1  Introduction 

With  ever  increasing  demands  placed  on  precision  and  re¬ 
liability  within  manufacturing  and  research  environments, 
magnetic  levitation  systems  are  finding  increased  utiliza¬ 
tion  within  such  applications  as  machine  tooling  due  to 
their  non-contact  (i.e.,  low  friction)  force  exertion  and  their 
capability  for  active  attenuation  of  mechanical  vibrations. 
However,  magnetic  levitation  systems  suffer  from  various 
control  complexities  such  as:  i.)  open-loop  instability,  ii.) 
inherent  nonlinearities  within  the  system  model,  iii.)  a  uni¬ 
directional  force  input,  and  iv.)  continuous  biasing.  In  re¬ 
sponse  to  these  difficulties,  researchers  have  employed  var¬ 
ious  modelling  and  control  development  techniques  to  ad¬ 
dress  the  aspects  of  i.)  through  iii.).  For  example,  Charara 
et  al.  [2]  constructed  a  nonlinear  model  of  an  inertial  wheel 
supported  by  active  magnetic  bearings,  for  which  a  sliding 
mode  controller  was  then  designed  to  stabilize  the  system. 
In  [7],  Levine  et  al.  proposed  a  nonlinear  feedback  con¬ 
trol  law  for  the  positioning  of  a  shaft  based  on  the  current 
complementarity  or  current  almost  complementarity  con¬ 
dition.  Queiroz  et  al.  [4]  utilized  a  nonlinear  model  of  a 
planar  rotor  disk,  active  magnetic  bearing  system  to  de¬ 
velop  a  global  exponential  position  tracking  controller  for 
the  full-order  electromechanical  system.  In  [8],  Mohamed  et 
al.  demonstrated  that  the  Q-parameterization  theory  could 
be  utilized  to  autobalance  the  rotor  of  a  vertical  shaft  active 
magnetic  bearing  system. 

Recently,  much  effort  has  been  directed  toward  the  area 
of  disturbance  suppression  within  magnetic  levitation  sys- 
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terns.  For  example,  in  [3],  Costic  et  al.  introduced  a  learn¬ 
ing  based  controller  to  asymptotically  regulate  a  magnetic 
bearing  system  while  compensating  for  periodic,  exogenous 
disturbances.  Rodrigues  et  al.  [9]  proposed  an  interconnec¬ 
tion  and  damping  assignment  Passivity-Based  controller  to 
yield  a  smooth  stabilizing  controller  for  the  active  magnetic 
bearing  system.  In  [6],  Gentili  utilized  a  model-based  regu¬ 
lation  approach  to  achieve  set  point  regulation  of  the  target 
with  the  disturbance  input  being  modeled  as  the  sum  of  a 
finite  number  of  sinusoids.  Similarly,  Behai  et  al.  [1]  de¬ 
signed  a  set  of  linear,  bounded-input  bounded-output  filters 
to  facilitate  the  utilization  of  standard  adaptive  techniques 
that  compensated  for  an  unknown  sinusoidal  disturbance 
signal.  Though  not  directly  targeted  at  magnetic  levita¬ 
tion  systems,  Xian  et  al.  [11]  proposed  an  adaptive  dis¬ 
turbance  rejection  approach  for  single-input  single-output, 
linear  time  invariant,  uncertain  systems  subjected  to  sinu¬ 
soidal  disturbances  with  unknown  amplitude  and  frequency. 
The  approach  of  [11]  utilizes  a  state  estimate  observer  in  a 
back  stepping  fashion  with  only  output  measurements  to 
achieve  asymptotic  disturbance  rejection. 

In  this  paper,  the  topic  of  disturbance  rejection  within  the 
magnetic  levitation  area  is  furthered  pursued  by  designing 
a  saturated  force  control  input  that  achieves  asymptotic 
target  position  regulation  despite  the  presence  of  a  nonlin¬ 
ear,  bounded,  periodic  disturbance.  The  control  develop¬ 
ment  differs  from  previous  disturbance  rejection  controllers 
in  that  restrictions  on  the  explicit  structure  of  the  distur¬ 
bance  signal  are  not  required  (i.e.,  the  disturbance  force 
need  only  be  bounded  and  periodic).  As  with  previously 
designed  magnetic  levitation  controllers,  the  proposed  con¬ 
trol  structure  must  contend  with  the  constraint  that  the 
actuation  force  is  unidirectional.  That  is,  the  magnetic  ac¬ 
tuator  can  only  exert  an  attractive  force  on  the  target  mass 
(the  earth’s  gravitational  field  is  utilized  to  produce  the  re¬ 
pulsive  action). 

The  remainder  of  the  paper  is  organized  as  follows.  The 
model  of  a  magnetic  levitation  system  actuating  on  a  tar¬ 
get  mass  suspended  in  the  gravitation  field  is  presented  in 
Section  II.  Section  III  identifies  the  control  objectives  and 
constraints  of  the  control  development.  The  saturated  con¬ 
trol  force  input  and  the  learning  based  disturbance  esti¬ 
mator  are  presented  in  Section  IV.  A  Lyapunov  stability 
analysis  is  utilized  to  illustrate  the  asymptotic  regulation 
of  the  target  position.  Simulation  results  are  presented  in 
Section  V. 
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2  System  Model 


u 


\  magnetic  levitation  system  consisting  of  a  target  mass 
suspended  vertically  in  the  gravity  field  subjected  to  a  non¬ 
linear  periodic  disturbance  force  can  be  modelled  by  the 
following  dynamics  [6] 

x  =  —  g  +  u  +  S  (t)  (1) 

where  x(t),  x  (/,),  x  (()  £  5ft1  represent  the  target  mass 
position,  velocity,  and  acceleration  signals,  respectively,  g 
—  9.81  m/s 2  £  3?1  denotes  the  gravitational  acceleration 
constant,  u 2  (t)  6  5ft1  represents  the  control  force  input1 , 
and  6  ( t )  €  Si1  denotes  the  nonlinear,  periodic  disturbance 
force  (note  that  the  system  of  (1)  has  been  normalized  with 
respect  to  the  target  mass  m). 


\J g  —  ki  tanh  (k^r)  —  6  ( t )  (5) 

where  8  (f)  £  5ft1  represents  a  learning  based  estimate  for 
8  ( t )  that  is  generated  on-line  via  the  following  expression 

S  ( t )  =  sat&0  (8  ( t  —  T))  +  &2  tanh  (fer)  (6) 

where  the  scalar  function  sats0  (•)  is  defined  in  the  following 
manner 

,  ,  .  f  e  for  lei  <  80  ,  . 

Sa  S°  ^  |  sgn  (e)  80  for  |e|  >  8a  ^ 

with  e  £  9ft1  representing  an  arbitrary  scalar  argu¬ 
ment,  sgn  (•)  denoting  the  standard  signum  function,  and 
ki,  k, 2,  k-3  £  9ft1  denoting  positive  scalar  control  gains  with 
fci  and  k2  selected  in  the  following  manner 

ki  +  fc2  <  g  -  80  (8) 


Remark  1  The  disturbance  force  8  (t)  is  assumed  to  be  pe¬ 
riodic  and  bounded  as  given  by  the  following 

8(t)  =  8{t-T),  \S(t)\<80<g,  (2) 

where  T  £  5ft1  denotes  the  known  period  and  60  represents 
a  positive  bounding  constant. 

3  Problem  Statement 

The  control  objective  is  to  design  a  suitable  force  input  that 
regulates  the  target  mass  position  to  a  desired  set  point  in 
the  presence  of  an  unknown,  bounded  periodic  disturbance. 
The  controller  design  is  complicated  by  the  presence  of  a 
nonlinear,  periodic  disturbance  and  also  by  the  fact  that  the 
levitation  system  can  only  exert  an  attraction  force  on  the 
target  mass.  In  order  to  facilitate  the  control  development, 
the  target  mass  tracking  error  signal  e  ( t )  £  5ft1  and  the 
filtered  tracking  error  signal  r  (t)  £  5ft L  are  definged  in  the 
following  manner 

e  =  x  —  Xd  r  =  e  +  ae  (3) 

where  Xd  £  9ft1  denotes  the  constant  desired  set  point  po¬ 
sition  and  a  £  5ft1  represents  a  positive,  constant  scalar 
control  gain.  The  proposed  control  force  is  developed  un¬ 
der  the  assumption  that  the  target’s  position  and  velocity 
signals  are  available  for  measurement. 

4  Saturated  Controller  Design 


Remark  2  The  structure  of  (6)  and  the  definition  of  (7) 
provide  a  bounded  disturbance  estimation  signal  8  (t)  in  the 
sense  that 

|^(f)|  <S0  +  k2; 

the  saturated  feedback  term  ki  tanh  (k^r)  in  the  control  input 
of  (5)  and  the  selection  of  the  control  gains  in  (8)  ensure 
that  the  radicand  of  (5)  is  non-negative. 

Remark  3  Based  on  the  definition  of  the  function  sats0  (-) 
in  (7),  the  following  inequality  can  be  illustrated  (See  Ap¬ 
pendix  A  for  details ) 

[<5  (t)  -  satSo  8  (t)] 2  <  [i5  (t)  -  6  (t)] 2  .  (9) 

Remark  4  Digital  implementation  of  the  learning  based, 
disturbance  force  estimator  8  ( t )  for  slowly  periodic  distur¬ 
bance  forces  (i.e.,  the  periodic  interval  T  is  large)  may  re¬ 
quire  significant  memory  allocation  if  the  control  sampling 
frequency  is  large.  That  is,  the  estimation  algorithm  of  (6) 
requires  the  logging  of  exactly  one  period  of  the  previous 
estimator  values.  For  example  if  the  controller’s  sampling 
frequency  is  adjusted  to  10  KHz  and  the  disturbance  force 
exhibits  a  period  ofT  =  1.0  sec,  then  the  selected  processing 
system  would  be  required  to  store  ( T )  ■  104  =  10,000  data 
points  in  memory.  Though  not  a  significant  concern  for 
memory  abundant  PC  applications,  this  phenomenon  more 
affects  the  compact,  memory  limited  DSP  or  microprocessor 
based  implementations  systems. 


The  development  of  the  regulation  controller  is  simplified 
by  rewriting  the  second  order  system  of  (1)  in  terms  of  the 
filtered  tracking  error  signal  r  ( t )  in  the  following  manner 

r  =  —  g  +  u2  +  8  (t)  —  a2 e  +  ar  (4) 

where  (1)  and  the  definition  of  (3)  have  been  utilized.  Based 
on  the  structure  of  the  ensuing  stability  analysis,  the  satu¬ 
rated  control  force  input  u 2  ( t )  is  designed  in  the  following 
manner 

1  Since  the  contribution  of  this  paper  lies  within  the  compen¬ 
sation  of  the  nonlinear  disturbance  force  8  (t),  we  have  chosen 
to  neglect  the  electrical  dynamics  for  the  magnetic  levitation 
system. 


After  substituting  the  control  force  input  of  (5)  into  the 
open  loop  dynamics  of  (4)  and  then  cancelling  common 
terms,  the  closed  loop  dynamics  for  r  (f)  can  be  formulated 
in  the  following  manner 

r  =  —ki  tanh  (far)  +  (<5  (t)  —  8  ( t ))  —  a2 e  +  ar.  (10) 


Theorem  1  If  the  control  gains  are  selected  to  satisfy  the 
following  condition, 


'  (2k\  +  fc2  —  Q3)  > 


|e(0)|+a|e(0)|  +  -e(0)2 


(ii) 
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where  e  (0)  and  e  (0)  denote  the  initial  target  posi¬ 
tion/velocity  tracking  error  signals,  then  the  control  force 
input  of  (5)  and  the  disturbance  estimator  6  (t)  of  (6)  en¬ 
sure  that  the  target  position  tracking  error  e  ( t )  is  driven 
asymptotically  to  zero  in  the  sense  that 


lim  e  ( t )  =  0. 

t — KX> 


(12) 


Proof:  In  order  to  illustrate  the  asymptotic  regulation  of 
e  (t)  in  (12),  the  following  non-negative  scalar  function  V  ( t ) 
is  defined  in  the  following  manner 


V  = 


Y  In  (cosh  (k3r))  +  i e 2 


t 

+  2 ^  J  [£(<r)-saf«0  (6  (cr))]2  da. 


(13) 


After  taking  the  time  derivative  of  (13),  substituting  in  the 

closed  loop  expression  given  by  (10)  and  utilizing  the  def¬ 
inition  of  (3),  the  following  expression  is  obtained  for  the 
time  derivative  of  V  ( t ) 

V  =  tanh  (k3r)  [— ki  tanh  ( k3r )  +  (6  ( t )  —  6  (f))] 

+  (ar  —  a2e)  tanh  ( k3r )  +  e  (r  —  ae) 

~wAs{t~T)~  sat6°  (* (t  -  T»]2  (14) 

+  2 ^  [S (*)  -  satgo  (Ht))]2. 

After  expanding  and  simplifying  the  terms  located  on  the 
second  row  of  (14),  V  ( t )  can  be  rewritten  in  the  following 
form 

V  =  tanh  ( k3r )  [—  k%  tanh  ( k3r )  +  (5  ( t )  —  6  (t))] 

+  (ar  —  a2e)  tanh  ( k3r )  +  e  (r  —  ae) 

+  2J-  |[<5(t)-sat«06(f)]2-  [(5(f)-6(t)]2j 
+2k2  tanh  (k3r)  (<5  ( t )  —  6(t )) 

-J-{kh^(k3r)}  (15) 

where  (2)  and  (6)  have  been  utilized.  The  inequality  of  (9) 
can  be  utilized  to  upper  bound  (15)  in  the  following  manner 


V  <  —  ^ki  +  tanh2  ( k3r ) 


+  tanh  ( k3r )  [  — a2e  +  ar]  +er  —  ae2. 


(16) 


In  order  to  complete  the  stability  analysis,  the  two  cases 
when  r(t)  =  0  and  r(t)  ^  0  are  examined. 

Case  1:  r(t)  =  0 

Based  on  the  fact  of  r(t)  =  0,  the  following  expression  is 
obtained 

tanh(fc3r)  =  0;  (17) 

thus,  the  time  derivative  of  V  ( t )  of  (16)  can  be  simplified 
into  the  following  form 


V  <  -ae2. 


(18) 


Case  2:  r(t)  ^  0 

Since  r  (t)  ^  0,  the  following  condition  exists 

tanh  (k3r)  /  0.  (19) 

Based  on  (19),  V  ( t )  of  (16)  can  be  rewritten  in  the  following 
manner 


i  ,  ,  fe  r 

*fcl+2  atanh(fer) 


tanh  ( k3r ) 


tann  ( k3r )  (20) 
e  •  tanh  (fc3r)] 


< 


~zT Bz  -  ^e2 


where  z  (t)  e  sIf2x  1  represents  an  auxiliary  vector  defined 
by  the  following 


and  the  matrix  B  ( t )  6 
manner 


;=  [  e  tanh  (for)  ]T  (21) 

,2x2  is  defined  in  the  following 


B  = 


a 

2 

a2 

T 


a 

~2 


2  tanh  (for) 


2  tanh  (for) 


,  fo 2  _  r 

1  2  tanh  ( k3r ) 


•  (22) 


In  order  to  ensure  that  the  time  derivative  of  V  ( t )  is  always 
negative  or  zero,  the  control  gains  must  be  selected  in  a 
fashion  such  that  the  matrix  B  of  (22)  is  positive  definite; 
thus,  the  control  gains  must  be  selected  to  guarantee  that 
the  following  conditions  are  satisfied 
a 


i-) 

ii.) 


>0, 


a  /,  ,  fo 
2  fcl  +  T 


tanh  (ksr) 
2 


(23) 


>  0. 


2  2  tanh  (£37-) 

Clearly  from  condition  i.),  the  control  gain  a  must  be  lower 
bounded  by  zero  in  the  sense  that  a  >  0.  Condition  ii.)  of 
(23)  can  be  expanded  and  simplified  to  obtain  the  following 


..  s  a  >.  a4  r2 

U  )  2(fcl  +  2^  4  4  tanh2  (k3r) 


>  0. 


After  utilizing  the  fact  that  [5] 

r{t) 


tanh  ( k3r ) 


<  | r  (f)|  + 


fo’ 


(24) 


(25) 


the  condition  of  (24)  can  be  rewritten  as  followings 

ii.)  a  (2 ki  +  k2-  a3)  >  ^|r|  +  j  .  (26) 

Based  on  the  structure  of  V(t)  defined  in  (13),  it  can  be 
seen  that 


V  >  —  In  (cosh  (for))  > 
k3 


H-Sh(2) 


(27) 


hence,  the  condition  of  (26)  can  then  be  rewritten  in  the 
following  form 

ii.)  a(2kx+k2-a3)  >  (V+^ln(2)  +  ^-)2.  (28) 
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With  selection  of  the  control  gains  according  to  i.)  of  (23) 
and  ii.)  of  (28),  the  auxiliary  matrix  B  is  positive  definite; 
thus,  V  (t)  of  (16)  can  be  upper  bounded  by  the  following 
expression 

V  <  —  ^e2  for  a  (2fci  +  fe2  —  a3)  > 


Prom  (29),  V(t)  <0  for  V/.  e  [0,  oo],  thus,  V (/.)  is  decreas¬ 
ing  or  constant  for  Vt  €  [0,  oo].  Therefore,  the  following 
sufficient  condition  for  (29)  can  be  obtained 

U<_|e2  for  a  (2ki  +  k2  -  a3)  > 

<3o) 

where  V (0)  is  explicitly  given  as  follows 

V(0)  =  t—  In  (cosh  (fcgr(O)))  +  ^e(O)2 

<  Ir(0)|  +  \e{Qf  {31) 

with  r  (0)  denoting  the  initial  filtered  tracking  error  signal. 
After  utilizing  the  inequality  of  (31),  the  expression  of  (30) 
can  be  rewritten  in  the  following  manner 
V  <  —  ^e2  for  a  (2ki  +  k2  —  a3)  >  (|e  (0)|  +  a  \e  (0)| 

+  Ie(0)2  +  ^ln(2)  +  ^)  (32) 

where  (3)  has  been  utilized. 

From  (13),  (27),  (32)  and  (18),  V  (t),  r(t),  e(t),  e(t)  6  £<*, 
and  e(t)  €  C2-  Based  on  the  previous  facts,  Barbalat’s 
Lemma  [10]  can  then  be  employed  to  illustrate  the  result 
presented  in  (12).  Standard  signal  chasing  arguments  can 
be  employed  to  illustrate  that  all  signals  remain  bounded 
during  closed-loop  operation. 

5  Simulation  Results 

In  order  to  illustrate  the  performance  of  the  saturated  con¬ 
troller  of  (5)  and  the  learning  based,  disturbance  estimator 
S (t)  of  (6),  the  system  dynamics  of  (1)  were  simulated  in 
Matlab’s  Simulink  simulation  package. The  desired  set-point 

position  was  selected  as 

xd  =  3  x  10“2  (m) ,  (33) 

and  the  periodic,  nonlinear  disturbance  signal  S  ( t )  is  shown 
in  Figure  2  (note  that  6  ( t )  is  bounded  and  exhibits  a  period 
of  T  =  1.0  seconds).  To  illustrate  the  convergence  abilities 
of  the  proposed  control  strategy,  the  target  was  set  to  the 
following  initial  position 

®(0)  =  -0.03  (m) .  (34) 

The  control  gains  were  adjusted  to  the  following  values  to 
obtain  the  best  performance 

fci  =  4  k2  =  1  k3  =  1  a  =  2  Sa  =  2.25 

(35) 

The  target  position  signal  x  ( t )  and  the  learning  based  es¬ 
timation  signal  S  ( t )  are  shown  in  Figures  1  and  2,  respec¬ 
tively. 


6  Conclusions 


In  this  paper,  a  saturated,  force  controller  targeted  for  mag¬ 
netic  levitation  systems  that  asymptotically  regulates  the 
position  of  a  target  mass  in  the  presence  of  a  nonlinear, 
bounded,  periodic  disturbance  has  been  presented.  The 
controller  development  does  not  require  restrictions  on  the 
explicit  structure  of  the  disturbance  other  than  an  upper 
bound  and  a  known  period.  In  addition,  the  magnetic  lev¬ 
itation  actuator  system  achieves  target  position  regulation 
despite  the  controller’s  limited  actuation  ability  in  that  only 
an  attractive  force  can  be  exerted  onto  the  target  mass  (i.e., 
the  controller  cannot  generate  a  repulsive  force  on  the  tar¬ 
get).  Future  research  in  the  area  of  disturbance  rejection 
for  magnetic  levitation  systems  will  focus  on  the  removing 
the  requirement  that  the  disturbance  period  to  be  known. 
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A  Inequality  Proof 

After  cancelling  the  resulting  matching  terms  from  the 
quadratic  expansion,  the  inequality  of  (9)  can  be  rewritten 
as  follows 

(sats08  (t))2  —  28  ( t )  ■  sats06  ( t )  <  (8  (t))”  —  26  ( t )  ■  8  (t) . 

'  -  2  -  i  (36) 

The  term  (6  (t)'j  —  28  ( t )  ■  6  (/,)  is  subtracted  from  both 

sides  of  (36),  and  the  resulting  expression  is  factored  to 
obtain  the  following  inequality 

(sats08  ( t )  —  8  (t))  (sat«05  ( t )  +8(t)~  26  (t))  <  0.  (37) 

In  order  to  verify  the  inequality  of  (37)  (and  thus  the  va¬ 
lidity  of  (9)),  the  remainder  of  the  proof  is  divided  into  the 
following  three  cases. 

A.l  Case  1:  |fi(t)|  <8„ 

Based  on  the  condition  that  |£(t)|  <  80,  the  function 
satsu6  ( t )  is  evaluated  as  follows 


sats06  (t)  =  8  (t) ;  (38) 

hence,  the  inequality  (37)  is  valid  for  the  conditions  of  Case 

1. 


A.2  Case  2:  6  ( t )  >  80 

Based  on  the  condition  that  8  (t)  >  So  and  from  the  defini¬ 
tion  of  (7),  the  function  sats06  (t)  is  evaluated  as 

sats0  6  (t)  =  So  (39) 

which  leads  to  the  following  two  conditions 

sats0S  (/,)  -  6  (t)  <  0,  satg06  ( t )  +  6  (t)  >  2 60  ;  (40) 

therefore,  the  inequality  of  (37)  is  validated  based  on  the 
fact  that  6  (t)  <  80. 

A.3  Case  3:  6  ( t )  <  —80 

In  a  similar  manner,  the  function  satg0  8  (<)  is  evaluated  as 
follows  for  the  condition  when  8  (t)  <  —So 

satSo6  ( t )  =  —  60  (41) 


which  leads  to  the  following  two  conditions 

sats06  (t)  —  8  (t)  >  0,  sats08  (f)  -I-  6  (L)  <  —26o 

(42) 

therefore,  the  inequality  of  (37)  is  valid  for  the  conditions 
of  Case  3  based  on  the  fact  that  <5  ( t )  >  — <5o  • 


Target  Position 


Figure  1:  Target  Position  x  ( t ) 


Learning  Based  Estimate  for  6{t) 


Figure  2:  6  ( t )  and  8  ( t ) 
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